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ABSTRACT Empirical studies of the first generation of stars and quasars will likely become 
feasible within the next decade in several different wavelength bands. Microwave anisotropy ex- 
periments, such as MAP or Planck, will set constraints on the ionization history of the intergalactic 
medium due to these sources. In the infrared, the Next Generation Space Telescope (NGST) will 
directly detect sub-galactic objects at Z J5 10. In the optical, data from the Hubble Deep Field 
(HDF) already places a constraint on the abundance of high-rcdshift quasars. However, the epoch 
of the first quasars might be first probed in X— ray bands, by instruments such as the Chandra 
X-ray Observatory (CXO) and the X-ray Multi-mirror Mission (XMM). In a 500 Ksec integration, 
CXO reaches a sensitivity of 2 X 10^^® erg Cm~^. Based on simple hierarchical CDM 
models, we find that at this flux threshold 10*^ quasars might be detectable from redshifts 
Z ^ 5, and ~ 1 quasar at Z ~ 10, in each 17' X 17' field. Measurement of the power spec- 
trum of the unresolved soft X-ray background will further constrain models of faint, high-redshift 
quasars. 
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1. INTRODUCTION 

The cosmic dark age ended when the first gas clouds condensed out of the primor- 
dial fluctuations at redshifts z = 10 — 20 (Peacock 1992; Rees 1996). These con- 
densations are hkely the sites where the first clusters of stars and the first quasar 
black holes appeared, giving birth to the first "mini-galaxies" or "mini-quasars" 
in the Universe. Despite the lack of observational data, this epoch has become a 
subject of intense theoretical study in the past few years. The recent interest can 
be attributed to forthcoming instruments: NGST could directly image sub-galactic 
objects at z ^ 10 in the infrared, while microwave satellites such as MAP or Planck 
could measure signatures from the reionization of the intergalactic medium (IGM). 
Currently, bright quasars are detected out to z ~ 5 (Fan et al. 1999). Although 
the abundance of optically and radio bright quasars declines at z ^ 2.5 (Schmidt 
et al. 1995; Shaver et al. 1996), a recent determination of the X-ray luminosity 
function (LF) of quasars from ROSAT data (Miyaji et al. 1998a) has not confirmed 
this decline. In this contribution, we point out that future X-ray observations might 
provide yet another probe of the first quasars and the end of the dark age at z ^ 10, 
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and that X-ray data might be uniquely useful in distinguishing quasars from stellar 
systems. 

2. THE APPEARANCE OF THE FIRST QUASARS AND STARS 

In popular Cold dark matter (CDM) cosmologies,the first baryonic objects appear 
near the Jeans mass (~ 10^ M©) at redshifts as high as ^; ~ 30 (Haiman & Loeb 
1999b, and references therein). At any rcdshift, the mass function of collapsed dark 
halos is given to within a factor of two by the Press-Schechter formalism. Following 
collapse, the gas in the first baryonic condensations is virialized by a strong shock 
(Bertschingcr 1985). Provided it is able to cool on a timcscalc shorter than the 
Hubble time, the shock-heated gas continues to contract. Depending on the details 
of the cooling and angular momentum transport, the gas then either fragments into 
stars, or forms a central black hole exhibiting quasar activity. Although the actual 
fragmentation process is likely to be rather complex, the average fraction /star of the 
collapsed gas converted into stars can be calibrated empirically so as to reproduce 
the average mctallicity observed in the Universe at z w 3. The observed ratio, 
inferred from CIV absorption lines in Lya forest clouds, is between 10"'^ and 10~^ 
of the solar value (Songaila 1997 and references therein). If the carbon produced 
in the early mini-galaxies is uniformly mixed with the rest of the baryons in the 
Universe, this implies /star ~2-20% for a Scalo stellar mass function. 

An even smaller fraction of the cooling gas might condense at the center of 
the potential well of each cloud and form a massive black hole, exhibiting quasar 
activity. In the simplest scenario, the peak luminosity of each black hole is propor- 
tional to its mass, and the light-curve, expressed in Eddington units, is a universal 
function of time. Indeed, for a sufficiently high fueling rate, quasars are likely to 
shine at their maximum possible luminosity, which is some constant fraction of the 
Eddington limit, for a time which is dictated by their final mass and radiative ef- 
ficiency. Here wc assume that the black hole mass fraction r = Mbh/Mhaio obeys 
a log-Gaussian probability distribution, p{r) = exp[— (logr — logro)^/2cr^], with 
logro = —3.5 and a = 0.5 (Haiman & Loeb 1999b). These values roughly reflect 
the distribution of black hole to bulge mass ratios found in a sample of 36 local 
galaxies (Magorrian et al. 1998) for a baryonic mass fraction of ~ (fib/^^o) ~ 0.1. 
We further postulate that each black hole emits a time-dependent bolometric lu- 
minosity in proportion to its mass, Lq = Mbh.fq = AfbhLEdd exp(— t/to), where 
i'Edd = 1-5 X 10^^ .Mbh/M© erg s~^ is the Eddington luminosity, t is the time 
elapsed since the formation of the black hole, and to = 10^ yr is the characteristic 
quasar lifetime. 

Finally, we assume that the shape of the emitted spectrum follows the mean 
spectrum of known quasars (Elvis et al. 1994) up to a photon energy of 10 keV. We 
extrapolate the spectrum up to ^ 50 keV, assuming a spectral slope of a=0 (or a 
photon index of -1). For reference. Figure 1 shows the adopted spectrum of quasars, 
assuming a black hole mass Mbh = lO^M©, placed at two different redshifts, Zg = 11 
and Zs = 6, and processed through the IGM, and assumed that reionization occurred 
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FIGURE 1. The observed spectra of quasars with a central black hole mass of Mbh = 10* Mq. 
The upper curves correspond to a source redshift of 2:3 = 6 and the lower curves to a source redshift 
of 2s = 11- In both cases we assume sudden reionization at Zr = 10. The dashed curves show the 
assumed intrinsic spectral shape, and the solid curves show the spectra after absorption by neutral 
intergalactic H and He. 



at z-c = 10 and that at higher redshifts the IGM was homogeneous and fully neutral. 
At lower redshifts, < 2; < Zr, we included the hydrogen opacity of the Lya forest 
given by Madau (1995), extrapolating his fitting formulae for the evolution of the 
number density of absorbers beyond z = 5 when necessary. As Figure 1 shows, 
the minimum black hole mass detectable by the ~ 2 x 10~^^ erg s~^ cm~^ flux 
Hmit of CXO (see below) is Mbh ~ 10* Mq at z = 10 and Mbh ^ 2 x 10^ M© 
at z = 5. In our model, the corresponding halo masses are Mhaio 3 x 10^^ Mq, 
and Afhaio 6 x 10^" M©, respectively. Although such massive halos are rare, their 
abundance is detectable in wide-field surveys. 



3. INFRARED: EXPECTED COUNTS WITH NGST 

The Next Generation Space Telescope]^ (NGST) will be able to detect the early 
population of mini-galaxies and mini-quasars. NGST is scheduled for launch in 
2008, and is expected to reach an imaging sensitivity of ^ 1 nJy (S/N=10 at spectral 
resolution A/AA = 3) for extended sources after several hours of integration in 
the wavelength range of 1-3.5/Ltm. Figure 2 shows the predicted number counts in 
the mini-galaxy and mini-quasar models described above, in a ACDM cosmology 
with (f]o,f^A,f^b,/i,cr8/i-i,") = (0.35, 0.65, 0.04, 0.65, 0.87, 0.96), normalized to a 

■'sec bttp://www. ngst.nasa.gov 
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FIGURE 2. Infrared Number Counts. The solid curves refer to quasars, while the long/short 
dashed curves correspond to galaxies with low/high normalization for the star formation efficiency. 
The curves labeled "5" or "10" show the cumulative number of objects with redshifts above z = 5 
or 10. 



5' X 5' field of view. This figure shows separately the number per logarithmic flux 
interval of all objects with redshifts z > 5 (thin lines), and z > 10 (thick lines). 
As the figure shows, NGST will be able to probe about ~ 100 quasars at z > 10, 
and ^ 200 quasars at z > 5 per field of view. The bright-end tail of the number 
counts approximately follows a power law, with dN/dF,^ oc F~'^-^ . The dashed 
lines show the corresponding number counts of mini-galaxies, assuming that each 
halo undergoes a starburst that converts a fraction of 2% (long-dashed) or 20% 
(short-dashed) of the gas into stars. These lines indicate that NGST would detect 
40 — 300 mini-galaxies at z > 10 per field of view, and ~ 600 — 10^ mini-galaxies 
at z > 5. Unlike quasars, galaxies could in principle be resolved if they extend over 
a scale comparable to the virial radius of their dark matter halos (Haiman & Loeb 
1997; Barkana and Loeb 1999). The supernovae and 7-ray bursts in these galaxies 
might outshine their hosts and may also be directly observable (Miralda-Escude & 
Rees 1997). Finally, we note that recent data in the J and H infrared bands from 
deep NICMOS observations of the HDF (Thompson et al. 1999) could already be 
useful to constrain mini-quasar and mini-galaxy models. 

4. OPTICAL: CONSTRAINTS FROM THE HUBBLE DEEP FIELD 

Although the infrared wavelengths are the best suited to detect the redshifted 
UV-emission from objects at z ^ 10, present data in the optical already yields 
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FIGURE 3. y-band number counts for high-redshift quasars. The lower panel shows differential 
counts for 3 < z < 3.2, as well as data adapted from Pei (1995). The upper panel shows cumulative 
number counts for z > 3.5. The top curve shows predictions of our standard ACDM model with 
Mjiij oc MjjjjIq. The bottom curve shows the counts in a model with halo circular velocities restricted 
to i)circ ^ 75 km s~^. This model is consistent with the lack of detections in the HDF. 

a constraint on quasar models of the type described above. The properties of faint 
extended sources found in the HDF (Madau et al. 1996) agree with detailed semi- 
analytic models of galaxy formation (Baugh et al. 1998). On the other hand, the 
HDF has revealed only a handful of faint unresolved sources, but none with the col- 
ors expected for high redshift quasars (Conti et al. 1999). The simplest mini-quasar 
model described above predicts the existence of ^ 10 B-band "dropouts" in the 
HDF, inconsistently with the lack of detection of such dropouts up to the ^ 50% 
completeness limit at y « 29 in the HDF. To reconcile the models with the data, a 
mechanism is needed for suppressing the formation of quasars in halos with circular 
velocities Vcirc ^ 50 — 75 km s~^ (see Figure 3 for the counts). This suppression 
naturally arises due to the photo-ionization heating of the intergalactic gas by the 
UV background after reionization (Thoul & Weinberg 1996; Navarro & Steinmetz 
1997). Alternative effects could help reduce the quasar number counts, such as a 
change in the background cosmology, a shift in the "big blue bump" component 
of the quasar spectrum to higher energies due to the lower black hole masses in 
mini-quasars, or a nonlinear black hole to halo mass relation; however, these effects 
are too small to account for the lack of detections in the HDF (Haiman, Madau & 
Loeb 1999). 
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FIGURE 4. The predicted X-ray luminosity function at 2 = 3.5 in our model (solid curves). The 
lower curve shows the effect of a cutoff in circular velocity for the host halos of fcirc > 50 km . 
The ROSAT data points are adopted from Miyaji et al. (1998a) and the dashed curve shows their 
fitting formula (for our background cosmology). 



5. X-RAYS: PREDICTIONS FOR CXO 

Quasars can be best distinguished from star forming galaxies at high redshifts by 
their X-ray emission. Detections of high-z quasars would therefore be highly valu- 
able: detections, or upper limits would help in answering the important question 
of whether the IGM at z ;^ 6 was reionized by stars or quasars, by yielding con- 
straints on the ionizing photon rate from high-z quasars. The simple quasar-model 
described above was constructed to accurately reproduce the evolution of the opti- 
cal luminosity function in the B-band (Pei 1995) at redshifts z ^ 2.2 (Haiman & 
Loeb 1998). However, it yields good agreement with the data on the X-ray LF, as 
demonstrated in Haiman & Loeb (1999c), and shown here in Figure 4. We regard 
this model as a minimal toy model which successfully reproduces the existing data, 
and use a straightforward extrapolation of this model to predict the X-ray number 
counts. In Figure 5, we show the predicted counts in the 0.4-6keV energy band of the 
CCD Imaging Spectrometer (ACIS) of CXO. Note that these curves are insensitive 
to our extrapolation of the template spectrum beyond 10 keV. The figure is normal- 
ized to the 17' X 17' field of view of the imaging chips. The solid curves show that 
of order a hundred quasars with z > 5 are expected per field at the CXO sensitivity 
of ~ 2 X 10~^^ erg s~^ cm~^ for a 5a detection of a point source. Note that CXO's 
arcsecond resolution will ease the separation of these point sources from background 
noise. The abundance of quasars at higher redshifts declines rapidly; however, a few 
objects per field are still detectable at z '-^ 8. The dashed lines show the results 
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FIGURE 5. The total number of quasars with redshift exceeding z = 5, z = 7, and 2 = 10 are 
shown as a function of observed X-ray flux in the CXO detection band. The soUd curves correspond 
to a cutoff in circular velocity for the host halos of J^circ ^ 50 km s~^, the dashed curves to a cutoff 
of fcirc ^ 100 km s~^. The vertical dashed line show the CXO sensitivity for a 5cr detection of a 
point source in an integration time of 5 X 10^ seconds. 

for a minimum circular velocity of the host halos of Wcirc > 100 km s""'^, and imply 
that the model predictions for the CXO satellite are not sensitive to such a change 
in the host velocity cutoff. This is because the halos shining at the CXO detection 
threshold are relatively massive, Mhaio ^ 10^^ Mq, and possess a circular velocity 
above the cutoff. In principle, the number of predicted sources would be lower if 
we had assumed a steeper spectral slope. For example, as figure 6 shows below, our 
model falls short of predicting the hard X-ray background, by about an order of 
magnitude at 10 keV. The difference could be explained by a change in our template 
spectrum to include a population of quasars with hard, but highly absorbed spectra 
(caused by the denser, and more gas rich hosts at high redshift). We note, however, 
that the agreement between the LF predicted by our model at z « 3.5 and that 
inferred from ROSAT observations would be upset by such a change, and require a 
modification of the model that would in turn tend to counter-balance the decrease 
in the predicted counts. 

6. THE X-RAY BACKGROUND 

Existing estimates of the X-ray background (XRB) provide another useful check 
on our quasar model. Figure 6 shows the predicted spectrum of the XRB in our 
model at z = (solid lines). The unresolved background flux is shown, obtained by 
summing the emission of all quasars whose individual observed flux at z = is below 
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FIGURE 6. Spectrum of the unresolved soft X— ray background in our model. We assume that 
the median X-ray spectrum of each source follows the mean spectrum of Elvis et al. (1994) up 
to 10 keV, and has a spectral slope of 0.5 (solid lines) or -0.5 (short-dashed lines) at higher 
photon energies. The lower curves show the spectra resulting from quasars with redshifts between 
2 < 2 < 5, and the upper curves include contributions from all redshifts z > 2. The long— dashed 
line shows the unresolved fraction (assumed to be 25%) of the soft X-ray background spectrum 
from Miyaji et al. (1998b). 

the ROSAT PSPC detection limit for discrete sources of 2 x 10~^^ erg cm~^ s~^ 
(Hasinger & Zamorani 1997). The short dashed hnes show the predicted fluxes 
assuming a steeper spectral slope beyond 10 keV (a = —0.5, or a photon index 
of -1.5). The long dashed line shows the 25% unresolved fraction of the soft XRB 
observed with ROSAT (Miyaji et al. 1998b; Fabian & Barcons 1992). This fraction 
represents the observational upper limit on the component of the soft XRB that 
could in principle arise from high-redshift quasars. As the figure shows, our quasar 
model predicts an unresolved flux just below this limit in the 0.5-3 keV range. The 
model also predicts that most ( ^ 90%) of this yet unresolved fraction arises from 
quasars beyond z = 5. The power spectrum of the unresolved background therefore 
might carry information on quasars at z > 5, and be useful in constraining the 
models (Haiman & Hui 1999, in preparation). The correlations in the background 
have recently been measured by Soltan et al. (1999, see also this Proceedings). 

7. DISCUSSION 

We have demonstrated that state-of-the-art X-ray observations could yield more 
stringent constraints on quasar models than currently available from the Hubble 
Deep Field (Haiman, Madau, & Loeb 1999). The X-ray data might provide the first 
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probe of the earliest quasars, complementing subsequent infrared and microwave ob- 
servations. More specifically, we have found that forthcoming X-ray observations 
with the CXO satellite might detect of order a hundred quasars per field of view in 
the redshift interval 5 ^ z ^ 10. Our numerical estimates are based on the simplest 
toy model for quasar formation in a hierarchical CDM cosmology, that satisfies 
all the current observational constraints on the optical and X-ray luminosity func- 
tions of quasars. Although a more detailed analysis is needed in order to assess the 
modeling uncertainties in our predictions, the importance of related observational 
programs with CXO is evident already from the present analysis. Other future in- 
struments, such as the HRC or the ACIS-S cameras on CXO, or the EPIC camera 
on XMM, which has a collective area 3-10 times larger than that of CXO, will also 
be useful in searching for high-redshift quasars. 

The relation between the black hole and halo masses may be more complicated 
than linear. With the introduction of additional free parameters, a non-linear (mass 
and redshift dependent) relation between the black-hole and halo masses can also 
lead to acceptable fits (Haehnelt et al. 1998) of the observed quasar LF near 2 ~ 3. 
Such fits, when extrapolated to higher redshift, can result in different predictions for 
the abundance of high-redshift quasars. From the point of view of selecting between 
these alternative models, even a non-detection by CXO would be invaluable. It is 
hoped further that either observations of the clustering properties of 2; ~ 3 quasars 
in the Sloan Digital Sky Survey, or a measurement of the power spectrum of the 
soft X ray background, would break model degeneracies (Haiman & Hui 1999, in 
preparation). 

Quasars emit a broad spectrum which extends into the UV and includes strong 
emission lines, such as Lya. For quasars near the CXO detection threshold, the fluxes 
at ^ 1/zm are expected to be relatively high, ~ 0.5-0.8 /iJy. Therefore, infrared 
spectroscopy of X ray selected quasars with the Space Infrared Telescope Facility 
(SIRTF) or NCST can identify the redshifts of the faint X-ray point-like sources 
detected by the CXO satellite. Such an approach could prove to be a particularly 
useful approach for unraveling the reionization history of the intcrgalactic medium 
at 2; <; 5. At present, the best constraints on hierarchical models of the formation 
and evolution of quasars originate from the Hubble Deep Field. However, HST 
observations are only sensitive to a limiting magnitude of y 29, and cannot 
probe the earliest quasars, beyond z ~ 6. The combination of X-ray data from the 
CXO satellite and infrared spectroscopy from SIRTF and NCST could potentially 
resolve one of the most important open questions about the thermal history of the 
Universe, namely whether the intergalactic medium was reionized by stars or by 
accreting black holes. 
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